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Abstract 

The Homestake result is about ~ 2a lower than the Ar-production rate, Qavi 
predicted by the LMA MSW solution of the solar neutrino problem. Also there is no 
apparent upturn of the energy spectrum (R = N Q b s /NssM) at low energies in SNO and 
Super-Kamiokande. Both these facts can be explained if alight, Amg x ~ (0.2- 2) -1CT 5 
eV 2 , sterile neutrino exists which mixes very weakly with active neutrinos: sin 2 2a ~ 
(10 -5 — 10 -3 ). We perform both the analytical and numerical study of the conversion 
effects in the system of two active neutrinos with the LMA parameters and one weakly 
mixed sterile neutrino. The presence of sterile neutrino leads to a dip in the survival 
probability in the intermediate energy range E = (0.5 — 5) MeV thus suppressing the 
Be, or/and pep, CNO as well as B electron neutrino fluxes. Apart from diminishing 
Qat it leads to decrease of the Ge-production rate and may lead to decrease of the 
BOREXINO signal and CC/NC ratio at SNO. Future studies of the solar neutrinos by 
SNO, SK, BOREXINO and KamLAND as well as by the new low energy experiments 
will allow us to check this possibility. We present a general analysis of modifications 
of the LMA energy profile due to mixing with new neutrino states. 



1 Introduction 



In the assumption of CPT invariance the first KamLAND result [1] and the results of SNO 
salt phase [2] confirm the large mixing angle (LMA) MSW solution of the solar neutrino 
problem [3, 4, 5]. Is the LMA solution complete? If there are observations which may 
indicate some deviation from LMA? 

According to the recent analysis, LMA MSW describes all the data very well [6, 7]: 
pulls of predictions from experimental results are below la for all but the Homestake ex- 
periment [7] . The generic prediction of LMA for the Ar production rate is 

Q Ar = 2.9-3.1 SNU, (1) 

which is about 2a higher than the Homestake result [8]. This pull can be 

• just a statistical fluctuation; 

• some systematics which may be related to the claimed long term time variations of 
the Homestake signal [8]; 

• a consequence of higher fluxes predicted by the Standard Solar Model (SSM) [9] 1 , 

• some physics beyond LMA. 

Another generic prediction of LMA is the "upturn" of the energy spectrum at low 
energies (the upturn of ratio of the observed and the SSM predicted numbers of events). 
According to LMA, the survival probability should increase with decrease of energy below 
(6 - 8) MeV [5]. For the best fit point the upturn can be as large as 10 - 15 % between 
8 and 5 MeV [10, 7]. Neither Super-Kamiokande (SK) [11] nor SNO [12] show the upturn, 
though the present sensitivity is not enough to make statistically significant statement. 

There are also claims that the solar neutrino data have time variations with small periods 
[13]. If true, this can not be explained in the context of LMA solution. 

Are these observations related? Do they indicate some new physics in the low energy part 
of the solar neutrino spectrum? In this paper we show that both the lower Ar-production 
rate and the absence of (or weaker) upturn of the spectrum can be explained by the effect 
of new (sterile) neutrino. The solar neutrino conversion in the non-trivial 3u- context 
(when the effect of third neutrino is not reduced to the averaged oscillations) have been 
considered in a number of publications before [4, 14]. In particular, modification of the v e - 
survival probability due to the mixing with sterile neutrino has been studied [15]. Here we 

1 For instance, the CWO-neutrino fluxes have rather large uncertainties. According to the SSM and in 
the LMA context they contribute to Qat about 0.25 SNU, so that reduction of the CNO- fluxes by factor 
of 2 (which is within 2cr of the estimated uncertainties) leads to reduction of the Ar-production rate by 
AQat ~ 0.12 SNU. 
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suggest specific parameters of the sterile neutrino which lead to appearance of a dip in the 
adiabatic edge of the survival probability "bath", at E — (0.5 — 2) MeV, and/or flattening 
of the spectrum distortion at higher energies (2 — 8) MeV. The dip suppresses the Be- (u e ) 
neutrino flux or/and other fluxes at the intermediate energies, and consequently, diminishes 
the Ar-production rate. It also diminishes or eliminates completely (depending on mixing 
angle and Am^) the upturn of spectrum. We comment on a possibility to induce time 
variations of signals by the presence of very small mixing with sterile neutrinos. 

The paper is organized as follows. In the Sec. 2 we introduce mixing with sterile neutrino 
and study in sec. 3 both analytically and numerically the conversion as well as the energy 
profile of the effect. In Sec. 4 physical consequences of the modification of the energy profile 
are considered. We calculate predictions for observables, the Ar-production rate, the Ge- 
production rate, the CC/NC ratio at SNO and the rate at BOREXINO, as functions of the 
mixing and mass of sterile neutrino in sec. 5. We consider an impact of the sterile neutrino 
on the global fit of the solar neutrino data in sec. 6, where we describe three possible 
scenario in sec. 6. In sec. 7 we discuss future checks of the suggested scenarios. We present 
a general analysis of possible modifications of the LMA profile by mixing with additional 
neutrino states in the Appendix. Our results are summarized in sec. 8. 



2 Sterile neutrino mixing and level crossing 

Let us consider the system of two active neutrinos, u e and u a , and one sterile neutrino, u s , 
which mix in the mass eigenstates u 2 and v$: 

vq = cos ct v s + sina:(cos# u e — sin# u a ), 
v 1 = cos a (cos 9 u e — sin6> u a ) — sin a u s , 

u 2 = sin 6 v e + cosQ v a . (2) 

The states v e and u a are characterized by the LMA oscillation parameters, 6 and Am 2 2 . 
They mix in the mass eigenstates v\ and v 2 with the eigenvalues mi, and m 2 . The sterile 
neutrino is mainly present in the mass eigenstate z/ (mass m ). It mixes weakly (sin a <C 1) 
with active neutrinos in the mass eigenstate v x 2 . We will assume first that m 2 > m > m x 
and consider the oscillation parameters of v s in the intervals: 

Am 2 ,! =ml-m\ = (0.2 - 2) x 10" 5 eV 2 , sin 2 2a ~ 10~ 5 - 3 • 10" 3 . (3) 



Let vi m , i>2mi u om be the eigenstates, and Ai, A2, Ao the corresponding eigenvalues of 
the 3z/-system in matter. We denote the ratio of mass squared differences as 

R, - (4) 

2 The introduction of mixing with the third active neutrino is straightforward. This mixing (if not zero) 
can produce small averaged oscillation effect and in what follows it will be neglected. 
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The level crossing scheme, that is, the dependence of \, (i = 0,1,2) on the distance 
inside the Sun (or on the density), is shown in fig. 1. It can be constructed analytically 
considering mixing of the sterile neutrino (the s-mixing) as a small perturbation. 

1). In the absence of s-mixing we have usual LMA system of two active neutrinos with 
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and \% MA has similar expression with plus sign in front of square root. Here V e = \f2Gp{n e — 
0.5n n ), and V a = — 0.5\/2GVn n are the matter potentials for the electron and non-electron 
active neutrinos respectively; n e and n n are the number densities of the electrons and neu- 
trons. For the sterile neutrino we have V s — 0. The 1-2 (LMA) resonance condition deter- 
mines the LMA resonance energy: 



E n 



Am^ cos 29 

We - V a ) " 



(6) 



2). Let us turn on the z/ s -mixing. In the assumption m\ < m < m 2 the sterile neutrino 
level \ s crosses \{ MA only. The level \2 MA essentially decouples. It is not affected by 
the s-mixing, and A 2 ~ \2 MA - Evolution of the corresponding eigenstate, z/ 2m , is strongly 
adiabatic. 



3). In general, the sterile level, X s , as the function of density, crosses \^ MA twice: above 
and below the 1-2 resonance density. Effects of the higher (in density) level crossing can be 
neglected since the neutrinos of relevant energies are produced at smaller densities. This 
can be seen in the fig. 1 where the second crossing of \{ MA and \ s would be on the left, 
if the density would continue to increase above the central solar density. Consequently, 
there are two relevant resonances in the system associated with 1-2 level crossing (the LMA 
resonance) and with 1-s crossing. For low energies (below the s-resonance) Ai ~ Af MA and 
Ao ~ A s . 

The Hamiltonian of the {y\^ A — v s ) sub-system can be obtained diagonalizing the u e — u a 
block of of the 2>v Hamiltonian, and then neglecting small 1-3 element. As a result 
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where Af MA is given in (5). The 1 — s resonance condition, 



^ MA (AmiJE,9,V e ,V a ) 
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determines the s-resonance energy 

^ 0.5m? + Am 2 ., cos 2 a 
E s = x — 3 x 



V e + V a 



x 1 ~ Ra ( 9 ) 

1 - 2R A + £cos2# + J (I - 2R A + £cos2#) 2 - AR A (l - i? A )(£ 2 - 1)' 



where £ = (V e — V a )/(V e + V a ) = n e /{n e — n n ). Notice that since Af MA is a non-linear 
function of the neutrino energy, the proportionality E s oc Am§j is broken and E s turns out 
to be complicated function of Am^. 

Another feature of the v\^ A — v s sub-system is that due to dependence of 9 rn on E, 
the effective mixing parameter in (7), oc sin2a:cos(# — 9 m ), also depends on the energy 
(decreases with E), though this dependence is weak. Indeed, in the case of small a and 
the s-resonance being substantially below the LMA resonance, we can take 9 ps 9 m in the 
first approximation, cos(# — 9 rn ) ps 1. Even in the LMA resonance, when 9 m = n/2, we get 
cos(# - 9 m ) = 0.97. 



3 Survival probability. Properties of the dip 



Let us find the v e survival probability. According to (2) the initial neutrino state can be 
written in terms of the matter eigenstates u im as 

v e = sm9^ n u 2m + cos6^(cosa^z/i m + sina^z/ 0m ), (10) 

where 9^ and oP m are the mixing angles in matter in the neutrino production point. 

Propagation of neutrinos from the production point to the surface of the Sun is described 
in the following way. v 2m evolves adiabatically, so that u 2m v i- Evolution of the two other 
eigenstates is, in general, non-adiabatic, so that 



A U V! + A 01 u , u 0m -»■ AxqUx + A QO v 0) 
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where A^ are the transition amplitudes which satisfy the following equalities: |v4 i | 2 — 
|A 10 | 2 = 1 — |A)o| 2 — 1 — l^n| 2 = P 2 . They can be found by solving the evolution equation 
with the Hamiltonian (7). P 2 is the two neutrino jump probability in the system v lm — v s . 

Using (10, 11) we can write the final neutrino state as 



v f = sin^z/ 2 e^ 2 + cos^ [cos a° m (A 11 u 1 + A 01 is ) + sina^A^ + A 00 u ) 
where <fi 2 is the phase acquired by v 2m . Then the survival probability is given by 



P ee = \(u e \u f }\ 2 w sin 2 9° m sin 2 9 + cos 2 9° m cos 2 9 [cos 2 oP m - P 2 cos 2a 



(12) 
(13) 
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Here we have neglected a small admixture of v e in vq\ (v e Wo) ~ 0. Also we have taken into 
account that the coherence of the mass eigenstates is destroyed on the way from the Sun to 
the Earth due to a spread of the wave packets and averaging effects. 

Similarly we obtain the transition probability of the electron to sterile neutrino: 

P es = \(v 8 \vf)\^cos 2 9l 



sin a + P2 cos 2a 



(14) 



Let us consider specific limits of the formula (13). If evolution is adiabatic in the s- 
resonance (which can be realized for the large enough s- mixing), we find P2 = and 

P ee = sin 2 9° m sin 2 9 + cos 2 9° m cos 2 9 cos 2 a° m . (15) 

In the opposite case of strongly non- adiabatic conversion (P 2 « 1) the probability equals 

P ee « sin 2 9° m sin 2 9 + cos 2 9° m cos 2 9 sin 2 a° m . (16) 

Notice that in spite of strong violation of adiabaticity in the s-resonance, the effect of 
s-mixing is still present due to the averaging of oscillations. 

The energy dependences of the probabilities can be easily understood using the results 
given in Eqs. (13 - 16). Let E a (n c ) and E s (n c ) be the LMA resonance energy and the s- 
resonance energy which correspond to the central density of the Sun n c . Then the following 
consideration holds. 

1) . For high energies, E > E a (n c ), neutrinos are produced far above the 1-2 resonance 
density, so that 9® m « ir/2. Then according to (13), P = sin 2 9, as in the 2v case, indepen- 
dently of properties of the s-resonance. The initial state coincides practically with z/ 2m , and 
the later propagates adiabatically. 

The s-resonance becomes operative at the energies of adiabatic edge, when 9® m deviates 
from 7r/2. This is the consequence of the fact that X s crosses the lowest LMA level Ai. 

2) . For low energies, E < E s (n c ), the s-resonance is not realized inside the Sun and 
s-mixing equals the vacuum mixing (cos 2 cP m ~ cos 2 a ~ 1). Then from (13) we get the 
usual adiabatic formula for the 2v case 

P e a f ab w sin 2 9° m sin 2 9 + cos 2 9° m cos 2 9 . (17) 



3). At the intermediate energies, crossing the s-resonance can be adiabatic (at E ~ 
E s (n c )), and moreover, the initial angle can be equal to ~ ir/2. Since the s-resonance is 
very narrow this equality is realized already at energies slightly above E s {n c ). In this case 
we get from (15) 

P ee w sin 2 9° m sin 2 9. (18) 
If also E < E a (n c ), so that 9° m w 6, the Eq. (18) leads to 

P ml n = Pee ~ Sin 4 9 . (19) 
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Pmin is the absolute minimum of the survival probability which can be achieved in the 
system. In general, P ee > sin 4 9, since E is not small in comparison with E a (sin 9® m > sin 9) 
and/or the adiabaticity is broken. 

For oP m m tt/2, which can be realized for E being slightly higher than E s , we find from 

(13) 

P ee pa sin 2 9° m sin 2 9 + cos 2 9° m cos 2 9 sin 2 a° m P 2 . (20) 

With the increase of energy the adiabaticity is violated, P 2 — > 1, and the probability ap- 
proaches the adiabatic one for the 2u system (17). 

In fig. 2 we show results of numerical computations of the v e survival probability P ee , 
and the survival probability of active neutrinos, (1 — P es ), as functions of energy. In our 
numerical calculations we have performed a complete integration of the evolution equations 
for the 3z/-system and also made averaging over the production region of the Sun. The 
analytical consideration allows us to understand immediately the numerical results shown 
in fig. 2. 

The effect of s-mixing is reduced to appearance of a dip in the LMA energy profile. A 
size of the dip equals: 

AP ee = P e L e MA -P ee = P es cos 2 9, (21) 

where P ee {E) LMA = P ee (E) adiah is the LMA probability given by the adiabatic formula (13). 
To obtain the last equality in (21) we used expressions for P ee from (13), P ee (E) LMA - from 
(17) and P es - from (14). Since cos 2 9 < 1 (the best fit value of LMA mixing, cos 2 9 = 0.714) 
according to (21) a change of the v e survival probability due to mixing with v s is weaker 
than the transition to sterile neutrino P es . The relation (21) is well reproduced in fig. 2. 

A position of the dip (its low energy edge) is given by the resonance energy taken at 
the central density of the Sun E s (n c ) (9). With increase of Am^ the dip shifts to higher 
energies. However, this shift is stronger than simple proportionality to Am^ as can be 
found from (9). For instance, the increase of R& from 0.1 to 0.2 leads to the shift of dip by 
factor 2.6 in the energy scale (see fig. 2). 

The maximal suppression in the dip depends on and a. For small R& (large spit 
between the two resonances) and large a (sin 2 2a > 10~ 3 ) the absolute minimum (19) can 
be achieved. Indeed, the condition for the minimum is nearly satisfied for the solid line in 
the upper panel of fig. 2 where P ee ~ 0.1. 

With increase of Ra (smaller split of the resonances) or/and decrease of a (stronger 
violation of the adiabaticity) a suppression in the dip weakens. Also with decrease of a the 
dip becomes narrower. 

Similarly one can consider crossing of A s with \i; MA . In this case the effect on P ee is 
weaker due to smaller admixture of v e in v 2 - Now the dip can appear at higher energies in 
the non-oscillatory part of the LMA profile where P ee ~ sin 2 9. 
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4 Observables and restrictions 



As follows from fig. 2, selecting appropriately the values of Ra and a (and therefore position 
and form of the dip) one can easily obtain significant suppression of Q Ar as well as the 
upturn of the spectrum (see fig. 3). There are, however, restrictions which follow from other 
experimental results. 

1). Ar-production rate versus Ge-production rate. A decrease of Q Ar is accompanied 
by decrease of Qc e (fig- 3). Since the LMA prediction for Q Ge is close to the central 
experimental value a possible decrease of Q Ge is restricted. Let us consider this correlation 
in details. 

The decrease of the Ar-production rate can be written as 

AQ Ar = Q Be r ■ AP e f + Qt • APjf + Q B r • AP e B e , (22) 

where Q% = 1.15 SNU, Qjgf = 0.64 SNU and Q B Ar = 5.76 SNU are the contributions to the 
Ar-production rate from the .Be-flux, the fluxes of the intermediate energies (pep, CNO) and 
the 5-neutrino flux according to SSM [9]. Here AP Be is the change of survival probability 
at E Be , AP™ f and AP B e are the changes of the effective (averaged over appropriate energy 
interval) survival probabilities for the intermediate energy fluxes and the boron neutrino 
flux respectively. 

The suppression of the Ge- production rate equals 

AQ Ge » Q B l ■ AP B : + QfSl ■ AP;f + Q B Ge ■ AP B , (23) 

where Q%\ = 34.2 SNU, ($£ = 11.7 SNU and Q% e = 12.1 SNU are the contributions to 
the Ge-production rate for the £?e-neutrino flux, the sum of pep- and CNO- fluxes, and the 
-B-neutrino flux correspondingly. AP Be is the same as in (22), whereas AF™' and AP B are 
approximately equal to those in (22). 

The changes of rates are correlated: 

AQ Ge = A(R A ,a)-AQ Ar , (24) 

where A is the constant which depends on the oscillation parameters. If the Be- (y e ) line 
is suppressed only, we would have A Be 30. If the neutrino fluxes at the intermediate 
energies are affected only, then A mt ~ 18, for the boron neutrino flux we find the smallest 
value A B ~ 2. 

In principle, the decrease of the Ge-production rate can be compensated by increase of 
the survival probability for the pp-neutrinos. This probability is given approximately by the 
average vacuum oscillations formula 

P ee (pp) w 1 -0.5 sin 2 20. (25) 
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From Eq. (25) it follows that the increase of P ee {pp) requires the decrease of mixing: 

Asin^ = APeM . (26) 

2^2P ee (pp) - 1 

The SSM contribution of the pp-neutrinos to Qc e equals Q P Q e = 69.7 SNU, therefore to 
compensate la (~ 5 SNU) decrease of Qae, one needs AP ee (pp) = 0.07. For this value of 
AP ee (pp) eq. (26) gives Asin 2 6> = —0.1. However, a decrease of sin 2 # is restricted by the 
high energy data (SK, SNO). Indeed, the survival probability for the boron neutrinos with 
E > 5 MeV is proportional to sin 2 9: 

P B ^asin 2 9, a « 1.1, (27) 

where the deviation of a from 1 is due to effects of the upturn and the v e regeneration in the 
matter of the Earth. So, the survival probabilities for the pp- and S-neutrinos are related: 

P ~ 1 _ 1 Pb + _|pj. (28) 
a a 1 

For the best fit value of mixing (sin 2 6> ~ 0.285) this equality gives AP PP « — 0.78APb. In 
turn, the survival probability Pb{> 5 MeV) is fixed by the CC/NC ratio: 



CC = P B 

NC 1-77 S (1-P B )' 



(29) 



where r] s is the sterile neutrino fraction in the state to which u e transforms. This relation 
does not depend on the original Boron neutrino flux. The solar neutrino data restrict 
i] s < 0.2, and therefore the presence of sterile component allows us to reduce the probability 
by a small amount only: AP ee ~ (0.2 — 0.3)i] s < 0.06. Moreover, according to fig. 2, the 
contribution of sterile neutrino to the high energy part of the spectrum is even smaller than 
0.2. 

2). The Ar-production rate versus the rates at SNO and SuperKamiokande. For large 
Ra and sin a the restriction appears from the charged current (CC) - event rate at SNO 
and well as from the rate of events at SK and the spectra. For free boron neutrino flux the 
suppression of Qa,- due to suppression of the boron electron neutrino flux can be written as 

AQ^ r = Qlf B AP e B e ' Ar , (30) 

where fs = Fb/F§ sm is the total boron neutrino flux in the units of the SSM flux. For the 
relative change of flux measured in CC-event, A[CC] = AF C c / Fq^ m we have 

A[CC] = f B AP e B e , (31) 

where AP B e is the change of the effective survival probability for the SNO energy range. 
With decrease of Qa t the rate [CC] decreases; we find 

A[CC] = 0.2Q Ar , (32) 
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and this relation does not depend of fs, so that for a given Qav, the decrease of [CC] can 
not be compensated by increase of 

Also the spectral information does not allow to strongly suppress Qat- 

5 Global Fit 

We have performed the global fit of the solar neutrinos data which takes into account the 
correlations of observables discussed in sec. 4. We use the same procedure of the fit as in 
our previous publications [7, 10]. In fig. 4 we show the dependence of the x 2 on Ra for fixed 
value of Am^. 

The following comments are in order. 

1. According to the fig. 4 the minimum Xo ~ 65.2 is achieved for 

R A = 0.9 sin 2 2a = 10~ 3 . (33) 

It corresponds to the unmodified Be- flux but suppressed pep- and CNO- fluxes. The upturn 
of the energy spectrum above 5 MeV is practically eliminated. Xo can be compared with 
X 2 ~ 66.6, for zero s-mixing. The improvements of the fit, A% 2 = 1.4, is not substantial. 
Notice however, that value Xo is n °t the absolute minimum. Furthermore, one should not 
expect significant improvement of the fit since the original pull was about 2a only, and 
quality of the global fit is very good in both cases. Finally with sterile neutrinos we have 
modified solution of solar neutrino problem with different set of predictions for observables. 

2. As follows from the fig. 4 certain regions of parameters of the sterile neutrino are 
strongly disfavored or excluded already by existing data. In particular, the region sin 2 2a = 
3 • 10~ 4 and Ra < 0.07 is excluded. It corresponds to strong suppression of the Be electron 
neutrino flux. 

In another strongly disfavored region: = 0.10 — 0.25, sin 2 2a > 10~ 3 , one has 

substantial suppression of the CC-signal at SNO and SK as well as distortion of the boron 
electron neutrino spectrum. For larger values of sterile neutrino mass, R& > 0.25, the dip 
shifts to higher energies and disappears. The conversion effects (and corresponding x 2 ) 
converge to the pure LMA solution case . 

6 Three scenarios 

Three phenomenologically different scenarios can be realized depending on the oscillation 
parameters, and therefore on the position and form of the dip. Three panels in the fig. 3, 
which correspond to different values of R&, illustrate these scenarios. Let us describe features 
of these three possibilities. 
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1) . Narrow dip at low energies: the _Be-line is in the dip. This corresponds to sin 2 2a < 
10" 4 and R A < 0.08 or 

0.5E Be < E s (n c ) < E Be , (34) 

where E Be = 0.86 MeV is the energy of the l?e-neutrinos (first panel in fig. 2 and solid line 
in fig. 3). The lower bound (34) implies that the pp-neutrino flux is not affected. In this 
case the .Be-line is suppressed most strongly; the u e fluxes of the intermediate energies (pep 
and CNO neutrinos) are suppressed weaker and the low energy part of the boron neutrino 
spectrum measured by SK and SNO is practically unaffected (see fig. 3). 

According to fig. 3 the value of coefficient in Eq. (24) A = 24. Taking the present la 
errors, 0.23 SNU and 5 SNU, for the Homestake and the combined Gallium result corre- 
spondingly, we find that the central experimental value of Qat can be reached at the price 
of the 2a decrease of Qce- 

The best compromise solution would correspond to sin 2 2a ~ 7 • 10~ 5 , when Qa t is la 
above the observation, and Q Ge is la below the observation. In this case the BOREXINO 
rate reduces from 0.61 down to 0.48 of the SSM rate (see sect. 7). 

For E s {n c ) being substantially smaller than E Be , the Be- line is on the non-adiabatic 
edge of the dip and its suppression is weaker. In this case larger values of sin a are allowed. 

As we have discussed in sec. 4 variations of the LMA parameters and the original boron 
neutrino flux do not allow us to compensate completely the changes of the observables 
(which worsen the fit) in the case when the Se-line is suppressed. 

2) . The dip at the intermediate energies: 

E Be < E s (n c ) < 1.4 MeV (35) 

(see the second panel in fig. 2 and the dashed lines in fig. 3). The .Be-line is out of the dip 
and therefore unaffected. A decrease of Qat occurs due to suppression of the v e components 
of the pep- and CNO- neutrino fluxes. 

In this case a decrease of Qa t is accompanying by smaller decrease of Q Ge in comparison 
with the previous case. For small enough mixing (so that the boron neutrinos are not affected 
strongly) we get from fig. 3 A = 15 in the relation (24). For larger sin 2 2a a suppression 
of the boron v e flux becomes substantial and A decreases further: A ~ 12. Now the value 
Qat = 2.8 SNU, which is la above the observation, can be achieved by just 0.4a reduction 
of Q Ge . 

The BOREXINO signal due to the Be- flux is unchanged, and also the observable part 
of the boron neutrino flux is affected very weakly. Change of the CC/NC ratio is about 
0.002. 

The optimal fit (see fig. 4) would correspond to sin 2 a = 10~ 3 , when Q Ar is diminished 
down to 2.75 SNU, at the same time Q Ge = 68 SNU and CC/NC = 3.22 in agreement with 
the latest data [2]. 
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3). The dip at high energies: 

E s (n c ) > 1.6 MeV (36) 

(see fig. 2, the panel for Ra = 0.2, and the dotted lines in fig. 3). Qat is diminished due 
to suppression of the low energy part of the boron neutrino spectrum. For sin 2 a = 1CT 3 . 
we find AQav = 0.17 SNU. At the same time a decrease of the Ge-pro duct ion rate is very 
small: AQc e ~ 0.5 SNU which corresponds to A = (2 — 3) in eq. (24). 

At sin 2 a = 10~ 3 there is already significant modification of the observable part of the 
boron neutrino spectrum and decrease of the total rate at SK and SNO. Also the CC/NC 
ratio decreases. According to fig. 3 at sin 2 2a = 10~ 3 , we have A(CC/NC) = 0.01. Fur- 
ther increase of Ra will shift the dip to higher energies, where the boron neutrino flux is 
larger. This, however, will not lead to further decrease of Qa t since the dip becomes smaller 
approaching the non-oscillatory region (see fig. 2). The BOREXINO signal (.Be-line) is un- 
changed. So, the main signature of this scenario is a strong suppression of the upturn and 
even a possibility to bend the spectrum down. 

Even for large Ra the influence of v s on the KamLAND results is negligible due to very 
small mixing. In contrast to the solar neutrinos, for the KamLAND experiment the matter 
effect on neutrino oscillations is very small and no enhancement of the s-mixing occurs. 
Therefore the effect of s-mixing on oscillation probability is smaller than sin 2 2a ~ 10~ 3 . 
For this reason the KamLAND result has not been included in the fit of data. 



7 Further tests 

How one can check the described scenarios? 

1) BOREXINO and KamLAND (solar) as well future low energy experiments [19, 20, 
21, 22, 23, 24, 25] can establish the suppression of the Se-neutrino flux in comparison with 
the LMA predictions, if the case 1) is realized. In BOREXINO and other experiments based 
on the z/e-scattering the ratio of the numbers of events with and without conversion can be 
written as 

RBorexino = -Pee(l r ) + f r Pesi (37) 

where r = a(u^e) / a(u e e) is the ratio of cross-sections. Using Eq. (21) we find an additional 
suppression of the BOREXINO rate in comparison with the pure LMA case: 

AR Borexino = Retina ~ R B ore X ino = (1 ~ r) AP ee + rP es W AP ee (l + T tan 2 6). (38) 

According to fig. 3, RB^exino can be diminished rather significantly. However, if the pre- 
diction for Qce is 2a (or less) below the experimental results, we find Rj£?rexino > 0-4 an d 
ARBorexino < 0.2. For the best fit value in the scenario 1): 

R L B™Lno~ 0.5, (ARBorexino ~0.1). (39) 
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Clearly, it will be difficult to establish such a difference. Furthermore, an additional sup- 
pression is mainly due to conversion to the sterile neutrino and the problem is to distinguish 
the conversion effect and lower original flux: the CC/NC ratio can not be used. Therefore 
not only high statistics results but also precise knowledge of the original fluxes is needed. 
The pep-flux is well known, however predictions of the CNO neutrino fluxes have larger 
uncertainties. 

2) . It may happen that the dip is at higher energies and the Be- flux is unaffected. In this 
case one expects significant suppression of the pep- and CNO- fluxes. Such a possibility can 
be checked using combination of measurements from different experiments which are sensi- 
tive to different parts of the solar neutrino spectrum. The radiochemical Li- experiment [26] 
has high sensitivity to the pep- and CNO- neutrino fluxes [27]. According to SSM [9], the 
C/VO-neutrino contribution to the -Be-production rate in this experiment is Q™ = 14.2 
SNU of the total rate Qgf ° = 52.3 SNU, so that Q C B ™° /Q Be = 0.27. For the CI- and 
Ga- experiments the corresponding ratios are substantially smaller: Q C j^° IQat — 0.05 and 
Q C G N e° IQce = 0.07. For the pep-neutrinos we get Q p ^/QBe = 0.176. 

Precise measurements of Q Be and Qc e and independent measurements of the B, pp and 
Be neutrino fluxes and subtraction of their contributions from Q Be and Qc e will allow to 
determine the CNO- electron neutrino fluxes. In general, to measure oscillation parameters 
and to determine the original solar neutrino fluxes one will need to perform a combined 
analysis of results from Ga-, CI-, Li- experiments as well as the dedicated low energy 
experiments [19, 20, 21, 22, 23, 24, 25]. Of course, new high statistics CZ-experiment would 
clarify the situation directly. 

3) . For i?A ~ 0.1 — 0.2 and sin 2 2a ~ 10~ 3 one expects significant suppression of the low 
energy part of the B- neutrino spectrum. As follows from figs. 5, at 5 MeV an additional 
suppression due to sterile neutrino can reach (10 - 15)% both in SK and SNO. The spectra 
with the s-mixing give slightly better fit to the data. Notice that there is no turn down of 
the SNO spectrum for R<\ = 0.2 and sin 2 #13 = 10~ 3 due to an additional contribution from 
the v — e scattering. Precision measurements of shape of the spectrum in the low energy 
part (E < 6 — 8 MeV) will give crucial checks of the described possibility. 

4) . In supernovae, neutrinos are produced at densities far above the LMA resonance 
density and propagation is adiabatic in the LMA resonance. So, even for very small 1-3 
mixing (sin 2 #13 > 10" 4 ) the adiabatic conversion u e — > z/ 2 is realized without any effect of 
sterile neutrino (as in the case described in Eq. (17)). If, however, the sterile level crosses 
the second level A 2 one may expect some manifestations of the sterile neutrino in the 
v e channel, provided that the mass hierarchy is inverted or the 1-3 mixing is very small 
(sin 2 # 13 < 10" 4 ). 

5) . Smallness of mixing of the sterile neutrino allows one to satisfy the nucleosynthesis 
bound: such a neutrino does not equilibrate in the Early Universe. For this reason sterile 
neutrinos also do not influence the large scale structures formation in the Universe. 
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6). A very small s-mixing means that the width of s- resonance is also very small. In 
the density scale An/n = tan 2a ~ 10~ 2 . Therefore 1% density perturbations can strongly 
affect conversion in the s-resonance [28] . If density perturbations (or density profile) change 
in time, this will induce time variations of neutrino signals. Since the effect of s-resonance 
is small, one may expect 10% (at most) variations of the Ga- and At- production rates. 

It seems that further precision measurements of the solar neutrino signals are the only 
possibility to check the suggested scenarios. 



8 Conclusions 

1. The low (with respect to the LMA prediction) value of the Ar-production rate measured 
in the Homestake experiment and/or suppressed upturn of the spectrum at low energies 
in SK and SNO can be explained by introduction of the sterile neutrino which mixes very 
weakly with the active neutrinos. 

2. The mixing of sterile neutrino leads to appearance of the dip in the survival probability 
in the interval of intermediate energies E = 0.5 — 5 MeV. The survival probability in the 
non-oscillatory and vacuum ranges is not modified (if sterile level crosses \{ MA ). 

3. Depending on value of i?A, that is, on a position of the dip, three phenomenologically 
different scenarios are possible: 

(i) the .Be-neutrino line in the dip; 

(ii) strong suppression of the pep- and CNO- neutrino fluxes, and the Se-neutrino line 
out of the dip; 

(iii) suppression of the boron flux only. 

The best global fit of the solar neutrino data corresponds to the unsuppressed fe- 
line, but strongly suppressed pep- and CNO- neutrino fluxes. Such a scenario requires 
sin 2 2a ~ 10~ 3 and Ra ~ 0.1. It predicts also an observable suppression of the upturn of 
the spectrum at SK and SNO. 

4. The present experimental results as well as relations between observables restrict sub- 
stantially possible effects of the dip induced by the s-mixing. 

5. The presence of s-mixing can be established by future precise measurements of the Be-, 
pep-, CNO- neutrino fluxes in BOREXINO [18] and KamLAND, as well as by measurements 
of the low energy part of the Boron neutrino spectrum ( < 5 — 6 MeV) in SNO. Study of 
the solar neutrinos seems to be the only possible way to test the scenarios described in this 
paper. There is no observable effects in laboratory experiment, as well as in astrophysics 
and cosmology. 

6. We have performed a general study of the effect of mixing with additional neutrino states 
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(see the Appendix). Only in the case when the sterile neutrino level crosses both the LMA 
levels, the effect of additional mixing can enhance the survival probability. This case is not 
realized, however, for additional sterile neutrino. In all other cases an additional mixing 
leads to suppression of survival probability. 

7. Even precise measurements of the high energy part of the solar neutrino spectrum may 
not be enough to reconstruct the energy profile of the effect at low energies. So, the low 
energy solar experiments are needed and they may lead to important discoveries. 

9 Note added 

Since the time we posted our paper on hep-ph, some new publications have appeared which 
are relevant for this study. 

1) . Lower value of the cross-section u N(p, 7) 15 measured by the LUNA experiment [29] 
leads to decrease of the predictions for the Ar-production rate are by AQ^ r = —0.1 SNU [30] 
(see our footnote 1). This reduces a difference of the LMA prediction and the Homestake 
result by about 0.5a. Notice that at the same time the Ge-production rate is dimished by 
AQ Ge = 2 SNU. 

2) . Larger values of the 7 Be(p,'-y) 8 B cross-section obtained in the recent measurements 
lead to significant increase of the predicted boron neutrino flux. Now the predicted flux is 
larger than than extracted from the NC event rate measured at SNO: f B = 0.88±0.04(ea;p)± 
0.23(theor) [31]. Being confirmed this may testify for partial conversion of the produced v e 
to sterile neutrino thus supporting scenario suggested in this paper. 
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Appendix: Profile of the effect and new neutrino states 

As we have established in the previous sections, mixing with sterile neutrino can modify 
the LMA energy profile, namely, suppress the survival probability in certain energy range. 
Here we present a general consideration of possible modifications of the LMA energy profile 
by mixing with new neutrino states. 

In general, an introduction of new states leads to decrease of P ee , since new channels 
open for disappearance of v e . What are conditions for increase of P ee ? 

In the LMA case the survival probability at low energies and at high energies are uniquely 
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related (28). So, in principle, measurements at high energies (E > 5 MeV) allow to recon- 
struct the profile at low energies provided that Am 2 is well determined. The latter can be 
achieved by KamLAND. Mixing with new states can change this high - low energy relation. 

In assumption that the coherence of all mass eigenstates is lost on the way to the Earth 
we can write the z/ e -survival probability after propagation in the Sun as 

Pee = J2 a ^\ 2 , (40) 

i 

where a, = |(^j|^/)| 2 is the probability to find the mass state % in the final state and 
U e i = (z/ e |z/j) is the mixing parameter. The quantities in eq. (40) satisfy the normalization 
conditions: 

5> = 1, ]T|f/ e ,| 2 = l. (41) 

i i 

At low energies, where neutrino conversion is due to the vacuum oscillations, the ad- 
mixtures of mass eigenstates are not changed and flavors are determined. Contributions 
from two different mass eigenstates add incoherently. In the 1v- case we have a\ = cos 2 9, 
a 2 = sin 2 9, U e i = cos 9, U e 2 = sin 9, and consequently, P ee = cos 4 9 + sin 4 9. 

The only way to increase P ee in vacuum, would be to restore the coherence (at least par- 
tially) of the two contributions, or decrease the mixing. In general, one should concentrate 
the electron flavor on one of the mass eigenstates and increase its admixture. 

Suppose additional neutrino states also produce the vacuum oscillation effect (no level 
crossing). Mixing of these new states with v e leads to decrease of |?7 e i| 2 or/and |t/ e2 | 2 , as 
well as di and a 2 , and one can easily show that 

P ee (2) > P ee (2 + n). (42) 

So, new states can lead to decrease of P ee only. 

Matter effects change a,. At high energies for the 2z/ case we get a\ ~ and a 2 ~ 1. Let 
\Uei\min and \U e i\ max be the largest and smallest mixing parameters correspondingly. Then 
it is easy to prove inequality 

\Uei\Ln < Pee < \U ei \ 2 max . (43) 

So, the only way to increase P ee is to change the admixtures of the mass states in such a 
way that aj, which corresponds to the largest \U e i\, increases. 

Let us consider one additional neutrino level (state) which mixes weakly with the LMA 
levels. Due to small mixing the LMA levels do not change significantly. If new (predom- 
inantly sterile) level crosses one of the LMA level only and P(i) is the probability that 
neutrino state does not transit to new level in this crossing, then 

P ee ^P(l) ai \U el \ 2 + a 2 \U e2 \ 2 . (44) 
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Here we put \U e o\ ~ 0. Since -P(l) < 1, the probability decreases as we have found in the 
specific case discussed in this paper. Similarly, if the new level crosses the second LMA 
level, the survival probability 

P ee ^a 1 \U el \ 2 + P{2)a 2 \U e2 \\ (45) 

decreases. Notice that if neutrino is produced far above the LMA resonance, so that a\ m 0, 
the probability equals P ee m P(2)a 2 |£/ e 2| 2 and for small P(2) (adiabaticity) the probability 
P ee can be strongly suppressed. 

To enhance P ee the new level should cross both LMA levels (in this case formulas above 
are not valid). If both crossings are adiabatic, the following transitions occur: 

V e « V 2m -> V 0m -> Vi m • (46) 

So that P ee = cos 2 6. If transitions are partially adiabatic, we find sin 2 6 < P ee < cos 2 9. 
Thus, the admixture of the mass state with the largest fraction of the electron neutrinos 
is enhanced. However, to get such a double crossing, the new level should have stronger 
dependence on the density than the dependence of the electron neutrino level. That is, the 
corresponding matter potential should be large: V x > V e . This is excluded: an additional 
sterile neutrino level can cross only one LMA level, thus leading to suppression of the survival 
probability. The results obtained in this paper are robust. 
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Figure 1: The level crossing scheme. The mass eigenvalues as functions of the distance from 

t 12 — 10 5 eV 2 and tan 2 9 = 0.4. The mass ratio is taken to 
Also shown is the position of 1-2 resonance (dashed vertical 
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E/Am 2 12 (MeV/eV 2 ) 

Figure 2: The survival probability of the electron neutrinos, P ee , (solid line) and survival 
probability of the active neutrinos, 1 — P es , (dashed line), as functions of E/Am\ 2 for 
different values of the sterile-active mixing parameter sin 2a. We take tan 2 9 = 0.4. Also 
shown is position of the 1-2 resonance for the central density of the Sun. (vertical dashed 
line). For Amj 2 = 7.1 • 10" 5 eV 2 the Be-\ine is at E/Amj 2 = 1.2 • 10 4 MeV/eV 2 , the pep- 
neutrino line is at E/ Am\ 2 = 2 • 10 4 MeV/eV 2 , the lowest (observable) energy, E = 5 MeV, 
and the highest energy of boron neutrino spectrum (~ 14 MeV) are at E/ Am\ 2 = 7 ■ 10 4 
MeV/eV 2 and 2 ■ 10 5 MeV/eV 2 correspondingly. 
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Am 21 2 =7.1x10" 5 eV 2 , tan 2 9=0.4 
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Figure 3: The Ar production rate (upper panel), the Ge production rate (second panel) the 
suppression factor for the BOREXINO signal and the CC/NC ratio at SNO as functions of 
sin 2 2a, for tan 2 9 = 0.4 and Am 2 x = 7.1 x 10~ 5 eV 2 . 
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Figure 4: The \ 2 °f the global fit of the solar neutrino data as a function of i?A for different 
values of the sterile-active mixing parameter sin 2 2a. We take tan 2 9 = 0.4 and Am 2 2 = 
7.1 x 10" 5 eV 2 , 
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Am 12 2 =7.1 10" 5 eV 2 , tan 2 0=O.4 and sin 2 2cc=10~ 3 
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Figure 5: The spectrum distortion (N osc /N SSM ) at Super-Kamiokande (upper panel) and 
SNO (lower panel) for different values of the mass ratio i?A and for the sterile-active mixing 
sin 2 2a = 10 -3 . The solid lines correspond to the pure LMA case (no sterile neutrino). 
Normalization of spectra have been chosen to minimize \ 2 fit of spectrum for each case. We 
show also the Super-Kamiokande and SNO experimental data points with statistical errors 
only. 
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